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It is known that any explicit averaging scheme of the type essential for describing the large scale 
behaviour of the Universe, must necessarily yield corrections to the Einstein equations applied in 
the Cosmological setting. The question of whether or not the resulting corrections to the Einstein 
equations are significant, is still a subject of debate, partly due to possible ambiguities in the averag- 
ing schemes available. In particular, it has been argued in the literature that the effects of averaging 
could be gauge artifacts. We apply the formalism of Zalaletdinov's Macroscopic Gravity (MG) 
which is a fully covariant and nonperturbative averaging scheme, in an attempt to construct gauge 
independent corrections to the standard Friedmann-Lemaftre- Robertson- Walker (FLRW) equations. 
We find that whereas one cannot escape the problem of dependence on one gauge choice - which 
is inherent in the assumption of large scale homogeneity and isotropy - it is however possible to 
construct spacetime scalar corrections to the standard FLRW equations. This partially removes the 
criticism concerning the corrections being gauge artifacts. For a particular initial choice of gauge 
which simplifies the formalism, we explicitly construct these scalars in terms of the underlying in- 
homogeneous geometry, and incidentally demonstrate that the formal structure of the corrections 
with this gauge choice is identical to that of analogous corrections derived by Buchert in the context 
of spatial averaging of scalars. 
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I. INTRODUCTION 

One of the basic assumptions of Cosmology is that 
the matter distribution in the Universe, when averaged 
on large enough scales, is homogeneous and isotropic. 
One then models the metric of spacetime on such large 
scales as having the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) form, i.e. - having homogeneous and 
isotropic spatial sections, and applies Einstein's General 
Relativistic field equations to determine the geometry of 
the Universe. The idea here is that just like the matter, 
the metric can also be "averaged" on large scales to yield 
a highly symmetric geometry. The averaging operation 
is usually an implicit (and often vague) assumption. It 
has long been known however that any explicit av- 
eraging scheme for the metric of spacetime and energy- 
momentum tensor of matter, must necessarily yield cor- 
rections to the Einstein equations - these equations be- 
ing ideally imposed on length scales where stars (and not 
galaxies) are pointlike objects. 

The problem of constructing an explicit and self- 
consistent averaging scheme for the Einstein equations in 
the Cosmological setting has a long history 0] ■ Recently, 
attention has been devoted to an averaging scheme de- 
veloped by Buchert [!, H| , who defines a spatial averaging 
operation for scalar quantities within a chosen 3 + 1 split- 
ting of spacetime, and derives effective equations which 
could have nontrivial effects on the assumed, homoge- 
neous and isotropic model of the Universe. The ques- 
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tion of whether or not the resulting corrections to the 
Einstein equations are significant, is still a subject of 
debate 0, |g, 0, @| and can be considered to be open. 
Buchert 's approach to averaging has also been debated 
on the grounds that it is concerned with spatial averages, 
and not spacetime averages. 

In such a context, it is perhaps desirable to work with 
an averaging scheme which starts from a spacetime aver- 
aging, and which allows us to construct the corrections 
to Einstein's equations as physically relevant, gauge inde- 
pendent objects that can (in principle at least) be obser- 
vationally tested. An averaging scheme for tensors on an 
arbitrary pseudo-Riemannian manifold has been devel- 
oped by Zalaletdinov and is detailed in Refs. [ij, fiol- fllj] . 
This scheme, called Macroscopic Gravity (MG), has been 
applied to Cosmology in some specific models [IH, [lj|, 
where the corrections were shown to be of the form of 
a spatial curvature term in the Friedmann equation. In 
this paper we address the following question : Can the 
equations of Zalaletdinov's MG, in a general setting and 
under suitable assumptions, be brought to a form where 
the corrections to the FLRW equations can be dealt with 
in a gauge independent manner? We will show that the 
issue of a 3 + 1 splitting can be reduced to one choice 
of gauge for the inhomogeneous metric which averages to 
the FLRW metric in its natural comoving spatial coordi- 
nates. While such a choice always needs to be made for 
any explicit computations (and is indeed inherent in the 
assumption of large scale homogeneity and isotropy), we 
will show that it is possible to at least formally construct 
spacetime scalar quantities which affect the FLRW equa- 
tions in a nontrivial manner. We will then explicitly con- 
struct these scalars assuming that the comoving FLRW 
spatial sections result on averaging in a particular gauge 
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choice which simplifies the formalism, and incidentally 
also demonstrate that for this gauge choice, the struc- 
ture of the corrections in Zalaletdinov's MG is formally 
identical to that of the corrections derived by Buchert [3| 
(see also Ref. 14]). Further, we will also show that in 
principle the allowed behaviour (with time) of these cor- 
rections is, in principle, more general than of the form of 
a spatial curvature, as was found in Refs. fl2l. [l3j. 

We have organized the paper as follows : For the ben- 
efit of readers unfamiliar with the details of Macroscopic 
Gravity, we have devoted considerable space to recalling 
this formalism in Section [TTJ which will also allow us to 
lay down our notation. Readers may wish to skim this 
section and go on to Section IIII| where we refine the 
assumptions of standard Cosmology and adapt them to 
the framework of MG, and also discuss the issue of gauge 
choice. Choosing a particularly convenient gauge for sim- 
plicity, we then explicitly compute the corrections to the 
FLRW equations in Section IIV1 and then we present the 
formal results for an arbitrary choice of gauge. This is 
followed by a comparison between the MG approach and 
the averaging scheme developed by Buchert, where we 
present the case that in principle the two approaches are 
expected to be identical in content, in the sense that 
they should yield the same physical results. We conclude 
in Section [V] by summarizing and discussing some addi- 
tional issues. 



II. THE MACROSCOPIC GRAVITY 
FORMALISM OF ZALALETDINOV 



sors on M. by boldface characters. For example, a. = 
(l/k\)ai 1 ...i k dx 11 A • • • A dx lk is a fc-form written in the 
coordinate basis of 1-forms {da; 1 }, while v = e a v a is a 
vector in the basis {e a }. The connection 1-forms are 
denoted by u> b i — (1%) dx- 7 where r£- are the Christof- 
fel symbols. For a fc-form p b .'.'. , we define the "exterior 
covariant derivative" D w associated with the connection 
uj a b , as follows : 

iX i>;; = dp a b ;:; - «* 6 a ,c: + u,\ a P f : + . . . . (1) 

The compatibility between the metric and the connection 
on M is expressed by the condition 

D^ffafc = dg ab - g a kU k b - gbk^ k a = > ( 2 ) 

with a similar condition for the inverse metric g ab . The 
Cartan structure equations are given by 

U3 a b A dx b = , (3a) 

du\ + u\Au c b = v\, (3b) 

where r a b is the curvature 2-form on M. which defines the 
Riemann curvature tensor via r a b = (l/2!)r a &cd dx c Ada; d . 
Finally, the structure equations ([3]) and the metric com- 
patibility condition ^ are supplemented by their re- 
spective integrability conditions, given by equations (I4al) . 
(lib)) and gc]) 

r a b A da; 6 = , (4a) 



Our aim in this section is to lay out the notation for 
the rest of the paper, and to recall the main results and 
concepts of Zalaletdinov's Macroscopic Gravity (hence- 
forth MG). The subject matter we will deal with is of a 
rather technical nature, and the results are in fact best 
presented in the language of differential forms. We will 
therefore begin with a brief recollection of some basic re- 
sults from differential geometry which will serve to fix 
the notation. We will then introduce additional nota- 
tion while describing the results of MG. The latter part 
of this section will essentially be a repetition of results 
obtained by Zalaletdinov Q and Mars and Zalaletdinov 
[l(i| . While we shall go through many of the algebraic de- 
tails in order to make this paper self-contained, we will 
not go into details of arguments justifying the various 
assumptions made in the averaging scheme, since our 
aim here is to apply the formalism developed in Refs. 

We shall deal with a differentiable manifold M. 
endowed with a metric g ab of Lorcntzian signature 
(— + ++)• We will use lower case Latin characters 
a,b,c, . . . i,j,k . . . to denote spacetime indices taking 
values 0,1,2,3; later we will also use upper case Latin 
characters A,B,C,... to denote spatial indices tak- 
ing values 1,2,3 in a chosen 3 + 1 splitting of space- 
time. We will denote all differential forms and ten- 



dr a b - u c b A r a c + cj a c A r c 6 = D w r a b = , (4b) 

9akr k b + g bk r k a = . (4c) 

We will now proceed to recall the main results of MG. 
The reader is referred to Refs. [^[Hl for details of deriva- 
tions of the formulae that follow. 

The construction of an averaged manifold M. from the 
manifold Ai begins with the definition of a bilocal exte- 
rior calculus on M, using a bilocal operator W" (x',x). 
Hereafter, the primed index refers to the point x' and 
the unprimed index to the point x, and the bivector 
Wj (x',x) transforms like a vector at x' and a coven- 
tor at x. This bivector is assumed to be idempotent [Iol |. 
i.e. W c a ,', {x\ x")Wf {x", x) = Wf {x', x); and to have the 
coincidence limit linV-^W" (x',x) = 6®. This ensures 
that Wj (x',x) has the inverse operator The 

bivector W° (x',x) serves three purposes - 
(a) It is used to "shift" the exterior derivative by eval- 
uating the derivative at x' and antisymmetrizing at x 
according to the rule (for say, a 2-form a) 

d' w a' = -a a /i,/ C '(x')Wi (x',x)dx j A dx a ' A dx b ' , (5) 
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where dx a is a coordinate basis 1-form at x', and so on. 
The shifted exterior derivative is used to define a fully 
bilocal exterior derivative d as 

d = d + d' w . (6) 

(b) W° is used to define the bilocal extension of tensors 
and &;-forms. For example, if P b (x) is a (1, 1) tensor on 
A4, then its bilocal extension is defined as 

P a b (x', x) = Wt (x', x)W a a , (x, x')Pi (x') . (7) 

Similarly, if a is a 2-form on M , its bilocal extension is 
defined as 

a(x',x) = -,a a >b>{.x')Wi (x\x)W h u(x',x)dx j Ad/. 
2! J 

(8) 

The definitions above can be easily generalised for an ar- 
bitrary tensor- valued fc-form on M. . Using the definition 
of bilocal extensions, the averages of say P b and a over 
a spacetime region X with a supporting point x, are then 
given by 

no*) = (n) ST = ^ J^x'^pux'^x) ■ 

a{x) = — { d 4 x' \/—g'a.{x' , x) ; Vs = [ d A x' \f—g' , 

(9) 

the subscript ST standing for 'spacetime'. (Note that 
a(x) is a (local) fc-form at x.) It can be shown 
that a necessary and sufficient condition for the bilocal 
exterior derivative to be nilpotent (dd =0) and hence 
for the average of a tensor-valued fc-form to be a single 
valued function of the supporting point x, is given by 

dW a ' = ; W a ' = wfdx 1 , (10) 

and further [l(| that a bivector W° satisfying all the re- 
quired properties always exists on an n-dimensional dif- 
ferentiable manifold with a volume n-form. 

(c) Finally, to completely define the exterior derivatives 
of averaged quantities, W°j is used to specify a Lie drag- 
ging of the averaging region X. This ensures that the 
volumes of the averaging regions constructed at nearby 
supporting points are coordinated in a well defined man- 
ner. Suppose x a and x a + £ a AA are the coordinates 
of two support points, where AA is a small change in 
the parameter along the integral curve of a given vector 
field £ a . Symbolically denoting the two points as x and 
x + £AA, the averaging region at x + £AA is defined in 
terms of the averaging region S at x, by transporting 
every point x' € X around x along the appropriate in- 
tegral curve of a new bilocal vector field S a defined as 
S a (x',x) — Wj (x',x)^(x), and thereby constructing 
the averaging region S(AA) with support point x + ^AA. 
(SeeRefs. [alS 

for further discussion on the significance 



of this averaging region coordination.) With this defini- 
tion, the exterior derivative of the average of a /c-form 
can be written as [9j 

dp£ = (dp£) ST + (div £ W A p£) ST - (div £ W) ST A p a b , 

where we have defined, in keeping with the notation of 
Ref. 1, 

div e W = Wf :a ,dx j 

= (W$ a , + Wfd a , In v^ 7 ) dx> . (12) 

Clearly, it is desirable to choose a coordination bivector 
Wj which satisfies 

div e W = 0, (13) 

since firstly, this allows us to commute the exterior 
derivative with the averaging according to 

dp a b = {dp a b ) ST , (14) 

and secondly, it implies that the volume of the averaging 
region is held constant during the coordination (To| , and 
is therefore a free parameter in the formalism. It can be 
shown [l0( that firstly, the general solution of Eqn. (fl7J|) 
for an idempotent coordination bivector is given by 

Wf(x',x) = ti(x')f- 1 J n (x), (15) 

where f^ l (x)d a — f m is any vector basis satisfying the 
commutation relations 

[f i ,f j ]=C k ij f k ■ C%= constant, (16) 

and secondly, that Eqn. (| 13|) with the coordination bivec- 
tor given by Eqn. (|15p is always integrable on a diffcr- 
entiable manifold with a given volume n-form. Further, 
for the special class of bivectors for which = 0, the 
vectors {f^} form a coordinate basis, with 'proper' coor- 
dinate functions 4> m (x) say, so that 

and satisfying Eqn. (| 13[) makes this proper coordinate 
system volume preserving, with g(4> m ) = constant. When 
expressed in terms of such a volume preserving coordi- 
nate (VPC) system, the coordination bivector takes its 
most simple form, namely 

Wf(x',x)\ piopeI =Sf. (18) 

Volume preserving coordinates in fact form a large class 
in themselves, generalizing the Cartesian coordinate sys- 
tem of Minkowski spacetime. For a discussion on the 
properties of VPCs and the associated bivectors W° , 
see Sec. 8 of Ref. [l(|- We now turn to describing the 
geometry on a manifold M. which is to be considered an 
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averaged version of the manifold M. . We denote the bilo- 
cal extension of the connection 1-form u) a b on M, as £l a b , 
so that 

Sl a b (x', x) = Tl, c , (x')W a a , (x, x')W b b ' {x',x)wi (x', x)dx c . 

(19) 

The key idea of MG is that the average Q a b of the con- 
nection 1-form on A4, is to be considered the connection 
1-form on the averaged manifold M. . The goal is then to 
average out the bilocal extensions of the structure equa- 
tions and the compatibility condition and their 
integrability conditions (j4]) , and to express them in terms 
of appropriate differential forms defined on M. . The bilo- 
cal extensions of Eqns. ([3]) and ([2]), are respectively given 
by 

n a b A dx b = , (20a) 
dfl a b + fl a c Afl c b ^r a b , (20b) 

Da gat = dgab - 9ak^ k b - gbk£l k a = , (20c) 

where, in the last equation, Dq is the bilocal covariant 
exterior derivative associated with the bilocal connection 
f2 a fc . The integrability conditions of Eqns. (|2H)l are given 
by the bilocal extensions of Eqns. Igj), 

r a b A dx b = , (21a) 
dr a b — Q c b A r a c + tt a c A r c b = , (21b) 

g ak r k b + g bk r k a = . (21c) 

To proceed with the averaging, a correlation 2-form is 
defined 

ZY, = (tt a fc A ^) ST - fl\ A . (22) 

The average of the curvature 2-form r a b on M. is denoted 
R a b = (r a b ) ST , and the curvature 2-form on the aver- 
aged manifold M. is denoted M a b , 

M a fc = dn a b + n a k a n k b . (23) 

Equations (|20a.|) and (|20b[) then average out to give 

fl a c A dx c = , (24a) 

M\ = K\ - Z\\ . (24b) 

It can further be shown [9j that the algebraic identities 
of the curvature 2-form on M. hold : 

M a c A da; c = ; M a a = , (25) 

To average out Eqn. (|21b|) . one needs to introduce a 
correlation 3-form which fixes the differential properties 
of the 2-form Z a fc l j- The differential properties of this 



3-form are then fixed by introducing a 4-form. However, 
accounting for these additional correlation forms is ex- 
tremely complicated. We will therefore make the simpler 
assumption (consistent with the formalism) of setting the 
correlation 3-form and 4-form to zero, with the condition 

DnZ» 6 V=0. (26) 

It will be an interesting exercise to check whether this 
assumption is, in fact, justified in case of averaging an 
exact inhomogeneous model. Eqn. (|2"6")1 has the integra- 
bility condition Q 

P (R a c A ZVj - ZVfc A R*,-) = . (27) 

Here the symbol P permutes the free indices in, 
say K a bj m n pairwise according to F(K a bj m n ) = 

(l/3!)(KV, m „ - KVy l „ + K im n \ -...), and any 
summed indices are ignored. Setting the correlation 
3-form and 4-form to zero also imposes the condition 

p(zV d AZ d / fc ) =0. (28) 

Eqn. (|21bp averages out to give the Bianchi identities for 
the curvature 2-form on M. 

D(jM [1 i) = dM a 4 - n k b A M° fc + n a k A M k b = . (29) 

To average out equations (I20c|) and (|21c|) . one needs to 
make additional assumptions. For a class of slowly vary- 
ing tensor fields (tensor-valued /c-forms) c™..' on M. such 
as the metric and other covariantly constant tensors, and 
Killing tensors, etc., if one assumes that 

(n\AcZ:) ST = fi\AcZ-: , (30a) 
{n a b a /\cZ:) ST - {n a b a ^) st ac™: , (30b) 

then Eqn. ()20c|) (and its analogue for g ah ) average out 
to give 

D s g ab = ; D^ a6 = 0. (31) 

Further, for a general slowly varying object c™..', the 
following identity holds 

(p 6 AC:') ST -R a 6 Ac:: 
-(n\A&ncZ\-) ST + n\AB n cZ: 

= -Z a h ™ A - . . . + ZY„ A c'jr + ■■■ , (32) 

which averages out Eqn. (|21c[) (and its analogue for g ab ) 
to give 

g ak M k b + g kb M k a = ; M\g kb + M a k g kb = . (33) 

Eqn. (|3"Tj) allows one to choose G ab = g a b, where G ab 
is the metric on the averaged manifold M.. In general 
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however, we have G ab ^ g ab , and one defines the tensor 
V ab = gab_ G ab to keep track Q f thig difference. (See Ref. 

Q for details.) However, we show in Appendix (|A lj) that 
when the averaged manifold is highly symmetric, as in 
the case of a manifold with homogeneous and isotropic 
spatial sections which we will consider, one finds that 
U ab = (see also the last paper in Ref. [Hi)- In the 
general case, it turns out that Eqn. (f3"2")l is all that is 
needed to average out the Einstein equations 

9 ak r kb - i*fa«r« = -< (mic) , (34) 

where k — 8ttGn, i?*-" 1 * ' is the microscopic energy mo- 
mentum tensor of the matter distribution, and the Ricci 
tensor r ab on A4 is defined according to the sign conven- 
tion r ab = T 3 gfoj ■ The averaging leads to the equations 

G ak M kb - l -8tG^M l3 = -k (t a b imic) ) ST 

+ (z a ijb - \s a b z% h )) r 

- (u ak M kb - i^'My) , 

(35) 

where M ab = M 3 ob . is the Ricci tensor on M. and we 
have defined 

Z a ij b — 2Z a ik k j b ; 7> a b j — Z a bm i j n dx m A dx n . (36) 

The averaged equations (f35j) differ from the usual Ein- 
stein equations by the correlation tensor which we define 

as 

CI = (z% b - ^% m ) g^-(u ak M kb - UtU^M^ 

(37) 

Hence, denoting the Einstein tensor on M. as E%, and 
defining the tensor T b via 

n = (t a b (mic) ) ST , (38) 

the averaged Einstein equations read 

El = -kTI + CI . (39) 

Since the left hand side of Eqn. (|39| is covariantly con- 
served by construction {E c b \ a = 0), where the semicolon 
denotes covariant differentiation with respect to the con- 
nection on M, in general one has 

(-/cZ? + Cg). o =0, (40) 

with no condition on and C^ separately. However, the 
condition (p?6")) after taking appropriate traces, reduces to 



(see the second of Ref. [TTj]) which implies that the av- 
eraged energy-momentum tensor is also covariantly 
conserved. 

It can also be shown that in 4 dimensions, the 720 a 
priori independent components of Z a bm l j n are subject to 
680 constraints arising from Eqns. (|27l) and ([28]) . This 
leaves 40 independent components which combine to give 
the 10 independent components of the correlation ten- 
sor Cg. The conditions in Eqns. (J27J) and (J28|) do not 
constrain the components of C£ , which follows from con- 
sidering the structure of those equations. (See also Ref. 

00 



III. A 3 + 1 SPACETIME SPLITTING AND THE 
SPATIAL AVERAGING LIMIT 

We are now in a position to apply the MG formalism 
to the problem of Cosmology. We start with the as- 
sumption that Einstein's equations are to be imposed on 
length scales where stars are pointlike objects (we denote 
such a scale as £mhom)- The averaging we perform will 
be directly at a length scale Lflrw larger than about 
100/i _1 Mpc or so. This averaging scale is assumed to 
satisfy i in hom < ^flrw < ^Hubble where i H ubbic is the 
length scale of the observable universe. The averaging 
will be assumed to yield a geometry which has homoge- 
neous and isotropic spatial sections. In other words, we 
will assume that the averaged manifold M admits a pre- 
ferred, hypersurface-orthogonal unit timelike vector field 
■D a , which defines 3-dimensional spacelike hypersurfaccs 
of constant curvature, and that v a is tangent to the tra- 
jectories of observers who see an isotropic Cosmic Back- 
ground Radiation. (These "observers" are defined in the 
averaged manifold - we will clarify below what they cor- 
respond to in the inhomogeneous manifold.) Throughout 
the rest of this paper, for simplicity, we will work with 
the special case where the spatial sections on M. defined 
by v a are flat. (In principle the entire calculation can be 
repeated for non-flat spatial sections as well.) One can 
then choose coordinates (t,x A ), A = 1,2,3, on M such 
that the spatial line element takes the form 

(M) ds 2 spatial =a 2 (t)S AB dx A dx B , (42) 

where Sab = 1 for 4 = B, and otherwise, and we 
have v a = (v 1 , 0,0,0) so that the spatial coordinates are 
comoving with the preferred observers. The vector field 
v a also defines a proper time (the cosmic time) r such 
that d T — v a d a — v*cV We will further assume that the 
averaged energy-momentum tensor Tg can be written in 
the form of a perfect fluid, as 

T a b = pv a v b +p^ bl (43) 

where the projection operator Tt b is defined as 



c a b;a = o, 



(41) 



< = 5 a b + v a v b , 



(44) 
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and p and p are the homogeneous energy density and 
pressure respectively, as measured by observers moving 
on trajectories (in M) with the tangent vector field v a , 

p = T a b v b v a ■ p=\ir b a T a b . (45) 

p and p are observationally relevant quantities, since all 
measurements of the matter energy density, especially 
those from studies of Large Scale Structure, interpret ob- 
servations in the context of the averaged geometry. (See 
Refs. [13, [3 for more careful treatments of this point.) 
An important consequence of the above assumptions is 
that the correlation tensor C b , when expressed in terms 
of the natural coordinates adapted to the spatial sections 
defined by the vector field v a , is spatially homogeneous. 
This is clear when the modified Einstein equations (|39|) 
are written in these natural coordinates. 

The existence of the vector field v a with the attendant 
assumptions described above, allows us to separate out 
the nontrivial components of the (FLRW) Einstein ten- 
sor E b on M. in a coordinate independent fashion - the 
Einstein tensor can be written as 

h(x) = E a b v b v a ; j 2 (x) ee 1 (TT b a EZ) , (46) 

where j\ [x) and ji ix) are scalar functions whose form de- 
pends upon the coordinates used. The remaining compo- 
nents given by n\E b v a and the traceless part of TT 1 a 'n h k E b , 
vanish identically. Since the energy-momentum tensor 
T b in Eqn. (|43[) also has an identical structure, this struc- 
ture is therefore also imposed on the correlation tensor 
C b . Namely, ^\C b v a and the traceless part of ^ l a ^\C b 
must vanish. This is a condition on the underlying in- 
homogeneous geometry, irrespective of the coordinates 
used on either M. or M, and is clearly a consequence of 
demanding that the averaged geometry have the symme- 
tries of the FLRW spacetime. (Alternatively, one could 
model the averaged energy-momentum tensor as having 
some anisotropic components, which would then have to 
be balanced by the corresponding components of the cor- 
relation tensor. We will not pursue this idea here, and 
will restrict ourselves to the more conservative assump- 
tion of a homogeneous and isotropic averaged matter dis- 
tribution P3"]) .) 

This leads us to the crucial question of the choice of 
gauge for the underlying geometry : namely, what choice 
of spatial sections for the inhomogeneous geometry, will 
lead to the spatial sections of the FLRW metric in the 
comoving coordinates defined in Eqn. ([4^)) ? Since the 
matter distribution at scale Linhom need not be pressure- 
free (or, indeed, even of the perfect fluid form), there 
is clearly no natural choice of gauge available, although 
locally, a synchronous reference frame can always be cho- 
sen. We note that there must be at least one choice of 
gauge in which the averaged metric has spatial sections 
in the form (|42p - this is simply a refinement of the Cos- 
mological Principle, and of the Weyl postulate, according 



to which the Universe is homogeneous and isotropic on 
large scales, and individual galaxies are considered as the 
"observers" travelling on trajectories with tangent v a . In 
the averaging approach, it makes more sense to replace 
"individual galaxies" with the averaging domains con- 
sidered as physically infinitesimal cells - the "points" of 
the averaged manifold M. This is physically reasonable 
since we know after all, that individual galaxies exhibit 
peculiar motions, undergo mergers and so on. This idea 
is also more in keeping with the notion that the Uni- 
verse is homogeneous and isotropic only on the largest 
scales, which are much larger than the scale of individual 
galaxies. (See also Ref. (l0( for a discussion of how this 
assumption of treating the averaging domains as being 
effectively point-like, is essential for the idempotency of 
the averaging operation.) 

Consider any 3 + 1 spacetime splitting in the form of 
a lapse function N(t,x J ), a shift vector N A (t,x J ), and 
a metric for the 3-geometry hAB{t,x J ), so that the line 
element on M. can be written as 

{M) ds 2 = — (N 2 — N A N A ) dt 2 +2N B dx B dt+h AB dx A dx B , 

(47) 

where Na = h,ABN B . At first sight, it might seem rea- 
sonable to leave the choice of gauge arbitrary. One could 
then formally consider a coordination bivector given by 
the Eqns. (|15[) and (fT?} , with x l denoting the coordinates 
in the chosen gauge and cj) m the VPCs; and demand for 
example, that the metric ((47|) (with say N A — 0) aver- 
age out to the FLRW form (with a nonsynchronous time 
coordinate in general). This would imply 

Goo = (.9oo) ST = ~/ 2 M ! G A = (goA) S T = ! 

Gab = (9ab)st = a 2 {t)8 A B ■ (48) 

There is no a priori reason to assume that the functions 
f(t) and a(t) are related. Note that the condition on 
the bilocal extension 170.4(2/, x) is in general nontrivial 
even when the components goA^x) are chosen to be zero. 
In the Appendix (|A 1[) we show that with the above as- 
sumptions, for a general lapse function N, the conditions 
T>Qg ab — (Eqn. (|3~lj) ) also allow us to choose 

U ij ee g ij - G ij = . (49) 

However, it turns out that for a general lapse function N, 
the explicit form of the correlation terms in the averaged 
equations with the assumptions (|48[) . is rather compli- 
cated (although one can write down these terms formally 
in a relatively compact manner). On the other hand, if 
we make the assumption that the spatial sections on M. 
leading to the spatial metric (|42]) on A4, are spatial sec- 
tions in a volume preserving gauge, then the correlation 
terms simplify greatly. This is not surprising since the 
MG formalism is nicely adapted to the choice of volume 
preserving coordinates. 

We therefore adopt the following procedure : We will 
first make the set of assumptions which allow us to use a 
volume preserving gauge. Using these assumptions in the 
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remainder of this section, we will introduce the notion of 
spatial averaging within the MG framework. In Section 
ITVlwe will calculate the correlation terms and display the 
modified equations resulting from this particular choice 
of gauge. This exercise can be thought of as a toy cal- 
culation with simplifying assumptions. Following that 
calculation we will show how the correlation terms can 
be generalized to the more physically relevant case where 
the gauge in the inhomogeneous metric is (formally) left 
unspecified. 

To begin our first calculation, we perform a coordinate 
transformation and shift to the gauge wherein the new 
lapse function N is given by N = 1 j\fh where h is the de- 
terminant of the new 3-metric 1~iab ■ In general, one will 
now be left with a non-zero shift vector N A ; however, 
the condition Nyh — 1 ensures that the coordinates we 
are now using are volume preserving, since the metric 



determinant is given by g 



-N h = — 1 = constant. 



We denote these volume preserving coordinates (VPCs) 
by (t, x) = (t, x A ) — it, x, y, z), and will assume that the 
spatial coordinates are non-compact. For simplicity, we 
make the added assumption that N A — in the inho- 
mogenous geometry [l5j ]. so that gn — —N 2 = — 1/h and 
9tA — 0. The line element for the inhomogenous manifold 
M becomes 



dP 



/i(t,x) 



+ h-AB(t, ~x)dx A dx 1 



(50) 



Note that in this gauge, the average takes on a particu- 
larly simple form : for a tensor pl(x), with a spacetime 
averaging domain given by the "cuboid" S defined by 

£ = {(f, x, y, z) | -T/2 < t < T/2, -L/2 <x,y,z< L/2} , 

where T and L are averaging time and length scales re- 
spectively, the average is given by 



{P)) ST (*' X ) = (Pi>ST ^' X ) 



J I ST v ' I Vj / ST 
^ r i+T/2 r +L/2 



TL 3 



dt' 



t-T/2 



dx'dy'dz' 



-L/2 



P i(t',x\y\z') 



(52) 



where the limits on the spatial integral are understood 
to hold for all three spatial coordinates. We define the 
"spatial averaging limit" as the limit T — > (or T <C 
^Hubble) which is interpreted as providing a definition of 
the average on a spatial domain corresponding to a "thin" 
time slice, the averaging operation now being given by 



Z3 



= / dx'dy'dz 1 



L/2 



p){t,x',y',z') 



0{TL-\ hUe ) 



(53) 



(Note the time dependence of the integrand.) Hence- 
forth, averaging will refer to spatial averaging, and will 
be denoted by (...), in contrast to the spacetime averag- 
ing considered thus far (denoted by (...) ST ). The choice 



of a cube with sides of length L as the spatial averag- 
ing domain was arbitrary, and is in fact not essential 
for any of the calculations to follow. In particular, all 
calculations can be performed with a spatial domain of 
arbitrary shape [16(. We will only use the cube for defi- 
niteness and simplicity in displaying equations. (An as- 
sumption essentially amounting to the limit in Eqn. (|53|) 
was also made in Ref. 13], in the context of averaging 
the spherically symmetric and inhomogeneous Lemaitre- 
Tolman-Bondi models.) The significance of introducing 
a spatial averaging in this manner is that the construc- 
tion of spatial averaging is not isolated from spacetime 
averaging, but is a special limiting case of the latter and 
is, in fact, still a fully covariant operation. 

For the volume preserving gauge, we now make the 
averaging assumption (with the averaging scale L = 
-^flrw) ? which reduces to 



G, 



(fJi 



Gab = (hAs) 



-1 



i 2 (i)S 



AB ■ 



(54) 

where a and / are some functions of the time coordi- 
nate alone. A few remarks are in order on this particu- 
lar choice of assumptions. Apart from the fact that the 
spacetime averaging operation takes on its simplest pos- 
sible form (f5"2")l in this gauge and allows a transparent 
definition of the spatial averaging limit, it can also be 
shown that the assumptions in Eqn. (|54|) are sufficient 
to establish the following relations : 



i 

jh) 



i 



(55) 



Here the second equality arises from the condition <f J = 
G u which can be assumed whenever the averaged met- 
ric is of the FLRW form (see Appendix (|A lj) ). The last 

equality follows on considering the conditions (j^bc) = 
(FLRW)pa^ j n ODV i ous notation, (the basic assumption of 
the MG averaging scheme) , details of which can be found 
in the Appendix (|A 2[) . Eqn. ([55")) reduces the line ele- 
ment on M. to the form 
dP 

d 2 (t)5 A Bdx A dx B . (56) 



^ds 2 



The line element in Eqn. (|56p clearly corresponds to the 
FLRW metric in a volume preserving gauge. In other 
words, the (spatial) average of the inhomogeneous geom- 
etry in the volume preserving gauge leads to a geometry 
with homogeneous and isotropic spatial sections, also in 
a volume preserving gauge. Note that the gauge in Eqn. 
(|5"6")) for the FLRW spacetime differs from the standard 
synchronous and comoving gauge, only by a redefinition 
of the time coordinate. The vector field v a introduced 
at the beginning of this section and which defines the 
FLRW spatial sections, is now given by 



,0,0,0) 



G ab V b 



-4,0,0,0, 



(57) 
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Before proceeding to the calculation of the correlation 
terms and the averaged Einstein equations, we briefly 
describe why it is important to consider the spatial av- 
eraging limit of the MG averaging operation. The key 
idea to emphasize is that an average of the homogeneous 
and isotropic FLRW geometry, should give back the same 
geometry. Since the FLRW geometry has a preferred 
set of spatial sections, it is important therefore to per- 
form the averaging over these sections. Further, since 
the FLRW metric adapted to its preferred spatial sec- 
tions depends on the time coordinate, it is also essential 
that the spacetime average should involve a time range 
that is short compared to the scale over which say the 
scale factor changes significantly. (See also Sec. 4 of Rcf. 
LL0|.) Clearly then, averaging the FLRW metric (denoted 
(FLRW) given in Eqn. (|56p (which is in volume pre- 
serving gauge) will strictly yield the same metric only in 
the limit T — ► 0. Namely, for the cuboid X defined in 
Eqn. HSU 

( {FLRW) 9a b ) = Hm ^ jf dfrfV (™><rf, x') 
= iFLRW) g ab , (58) 

which should be clear from the definition of the metric. 
The result (( FLRvr >g ab ) = {FLRW) g a b in the spatial aver- 
aging limit can also be shown to hold for the FLRW met- 
ric in synchronous gauge, where the coordination bivec- 
tor Wj can be easily computed using the transformation 
from the VPCs (i, x A ) to the synchronous coordinates 
(r,y A ) given by 

^S'W) • 9 " W ' <59) 

The transformation (|59|) will also later allow us to write 
the averaged equations in the synchronous gauge for the 
averaged geometry. 

We now proceed to calculating the correlation 2-form 
Z a b l J - and thereby the averaged Einstein equations. 

IV. THE CORRELATION 2-FORM AND THE 
AVERAGED FIELD EQUATIONS 



where the dot will always refer to a derivative with re- 
spect to the VPC time t, and h AB is the inverse of the 
3- metric hAB- (This also gives the symmetric tensor 
Qab — (l/2N)hAB, which is the negative of the extrinsic 
curvature tensor.) The traceless symmetric shear tensor 
<j£ and the shear scalar er 2 are defined by 



*i = ei-(Q/3)8i ; ^ = 1^, (61) 



where = 9^ = (1/N)dt \n^fh is the expansion scalar. 

The connection 1-forms u> l j = r* fe dx' c can be easily 
calculated in terms of the expansion tensor, for an arbi- 
trary lapse function N . Specializing to the volume pre- 
serving gauge (N — /i -1 / 2 ), the bilocal extensions of 
the connection 1-forms are trivial and are simply given 
by 



n i j (x',x)=r i jk (x')dx k . (62) 



Since G a b — g a b, the connection 1-forms O l j for the aver- 
aged manifold M. are constructed using the FLRW metric 
in volume preserving gauge given in Eqn. ([56]), and can 
also be easily evaluated. 

We can now construct the correlation 2-form Z a b * • de- 
fined in Eqn. (|22[) . For completeness, we will display all 
the nontrivial components Z a b ^ , although not all of them 
will be relevant for the final equations. The condition 
N = h^ 1 / 2 has the effect that several of the Christoffel 
symbols become related to each other. For example, we 
have rg = -dt(lnVh) = -T A A = -(1/Vh)@, and so on. 
In the following, ^r^ c denotes the Christoffel symbol 
built from the 3- metric Hae, and we have defined 



A. Results for the Volume Preserving Gauge 

We start by defining (in any gauge with N A = 0) the 
expansion tensor 0^ by 



a dt 



(63) 



e A 



B= 2N h hcB ' 



(60) 
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We have, 



7O _ 17 

^ A — A A 



(@@aj) + (d A (\nVh)dj{lnVh)j - 3a a H 2 S A j\ dx J A dt + (yh@ AJ d K (lnVh)) dx J A dx K , (64a) 



tO A 



tA 



~h AK d K (lnVh)dj{lnVh)J - 3H 2 Sj dx J A dt + (^=<djd K (hi Vh)J dx J A dx K , 



(64b) 



rvO A 

& B 



rjA 

^ B 
1 



_0(3)p^ 



-Le B dj(inVh) 



dx J A dt + ( P)r£j8 K Qj3. Vh)) dx J A d 



J< 



(64c) 



70 

^ A B 



2(Vhd [A (lnVh)Q B]J )] dx J A dt + [(/ie A/ e SK ) - a 16 # 2 <W<W] dx J Ad 



7 A A x K 



(64d) 



rrB 

A 



^=d A (]nVh)Q B 



-^=h BK d K (\nVh)0 Au 



dx J A dt + [(QajQk ) ~ a*H 2 8 AJ 5 B : ] dx J A da; K , 

(64e) 



rjQ B _ rjB 

^ A C* ~ — ^ C A 



9 A (ln vft) (3) rgj) + <6Ajeg) - a 8 H 2 8 AJ 8 B 1 ] dx J A di+ (VhQ AJ {3) T^ K ) dx J A dx K , (64f) 



T A B 
J 



dx J A dt 



7 Vj V K 



H 2 5j5 B K 



dx J A dx K . 



(64g) 



rA B 

J C 



7 B A 

& c 



I^^(in^) <3>rgA + / Ie^eg\ - H 2 5js B 



dx J A dt + <^^ e J (3, rcjfJ> dx J A da; A ' , (64h) 



rjA C 

A B D 



Vh 



ftC (3) r A 



J_ B A (3)pC 



dcc J A dt 



(3) rV 3) rg K W 7 Ad/ 



(64i) 



r 



where we have used the relation ^T BJ — deQnv^i), 
and the square brackets around indices indicate antisym- 
metrization (p [y] = (l/2!)(py ~Pji). ) 

It is now straightforward to use the relations in Eqn. 
(|36[) (note the unconventional normalization of the 2- 
form) to read off the components Z a bm l - n and hence per- 
form the required summations to construct Z a i - b . This, 
together with the relation g ab = G ab , allows us to con- 
struct the correlation tensor C b defined in Eqn. (pTT)) 



h — \ 6 ijb l^ b Zj ijm 



(65) 



Now, the components of the Einstein tensor for the 
averaged spacetime with metric (I56p are given by 



(66) 



E\ = 3a 6 H 2 ; E A = = Ef , 



TpA _ =6 rA 

E B - a d B 



2 \ - + 3H 

a 



2 \ ,. _g-2 



where the peculiar splitting of terms in the last equation 
is for later convenience. Recall that the overdot denotes 
a derivative with respect to the VPC time t, not syn- 
chronous time. In terms of the coordinate independent 
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objects introduced in Eqn. (|46|) , we have 



-3a e H 2 



32(x) - a 6 



2 ( - + W 

a 



(67) 

^From the averaged Einstein equations in ()39|) we next 
construct the scalar equations which in the standard case 
would correspond to the Friedmann equation and the 
Raychaudhuri equation. These correspond to the Ein- 
stein tensor components, 

E a b v b v a =j 1 (x) ■ ir b a E a b +E a b v b v a = 3j 2 (x)+j 1 (x), (68) 

and are given by 



3a 6 H 2 = {KTl~C a b )v a v b 
1 

= up — — 



(69a) 



6a° [ - + 3H' 

a 



(_ kT £ + C a b ) (v a v b 



-k (p + 3p) 



(69b) 



Here Eqn. (l69a|) is the modified Friedmann equation 
and Eqn. (|69bj) the modified Raychaudhuri equation (in 
the volume preserving gauge on .M). We have used Eqn. 
(|45p , with the overbar on p and p reminding us that they 
are expressed in terms of the nonsynchronous time t, and 
we have defined the correlation terms 



Q 



(1) - 



^)-> 2 ))-*(r 2 

1 



^( F 6 2 ) = (3Hf , (70a) 



rAB 



(3) r J (3) r C 
1 AC 1 B.J 



- (d A (lnVh)d B (lnVh))] , (70b) 



G w = 5 6 (±eiez) - ±s AB (e AJ e J B ) , ( 70c ) 



a 6 /~h AB d A {lnVh)d B (ln Vh)\ 



5 AB (d A (lnVh)d B (lnVh)) . (70d) 



5< 2 



In defining we have used the relation 2 — S^B^ = 
(2/3)6 2 - 2cr 2 . QW and Q (2) are correlations of the 
extrinsic curvature, whereas and are correla- 
tions restricted to the intrinsic 3-geometry of the spa- 
tial slices of M.. Since the components of C b are not 
explicitly constrained by Eqns. (|2"T)) and l[28p. we can 
treat the combinations (1/2)(QW + S {1) ) = -Cg and 



Cq) as independent, sub- 



2 (Q (l) + Q ( 2 ) +5 ( 2 ))-^ A -^ ; 

ject only to the differential constraints (|41[) which we will 
come to below. 

As discussed in the beginning of Section IIII1 the re- 
maining components of C b must be set to zero, giving 
constraints on the underlying inhomogeneous geometry. 
In coordinate independent language, these constraints 
read 



= <C a b v b 



o. 



(71) 



Eqns. (fTTj) reduce to the following for our specific choice 
of volume preserving coordinates, 



c° A = o 



c A = o 



c A B -\s A B (c J j) 



0. 



(72) 



It can be shown that the VPC assumption N — h 1 / 2 
reduces the correlations and defined in Eqns. 
(|70c|) and (|70dj) . as well as several terms in the explicit 
expansion of Eqn. (|72|) . to the form 



(.9oo ) 



AiJ-niii j-iii \ / AB\ /-pi' 

5005 r 6ici r ilfcl ^ - (g ){r b2 



A>) ■ (73) 



The assumption in Eqn. (|30b[) (which is fundamental to 
the MG formalism) now shows that one can write 

(g og AB n c r Jk ) = (goog AB ) (r a bc r Jk ) 

' h AB \ 

(n c r 3k ) ■ (74) 

An interesting point is that the VPC assumption 
N = h~ ll>2 further allows us to assume (h AB /h) = 
(h ) (1/h) consistently with the formalism (details in 
Appendix (|A 2|> ) . Using Eqn. ([55)1 this gives us 

(ir)-? <»">• <75 » 

This shows that the correlation terms and in 
fact vanish, 



Q( 2 > = = , 



(76) 



and leads to some remarkable cancellations in Eqns. (|72[) . 
which simplify to give 



(VhQj B 



(3)pB 
1 AK 



(VhO 



IK 



(3) r B AB ) 



5 J'<(±=q b ^t 



rAJ 



Vh K " 



.IB 



o, 

(77a) 

-0, 
(77b) 



6 JK {^V A C ^T C KB ) - 5 AJ 



(3) r C (3) r K 
1 JC 1 BK 



(77c) 
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These simplifications are solely a consequence of assum- 
ing that the inhomogeneous metric in the volume pre- 
serving gauge averages out to give the FLRW metric in 
standard form. In general, these simplifications will not 
occur when the standard FLRW metric arises from an 
arbitrary choice of gauge for the inhomogeneous metric. 

In order to come as close as possible to the standard ap- 
proach in Cosmology, we will now rewrite the scalar equa- 
tions (|69|) (which are the cosmologically relevant ones) 
after performing the transformation given in Eqn. (|59|) 
in order to get the FLRW metric to the form 

{M) ds 2 = -dr 2 + a 2 (T)5 A Bdy A dy B ; o(r) = a(t{r)) . 

(78) 

Since Eqns. (|69p are scalar equations, this transforma- 
tion only has the effect of reexpressing all the terms as 
functions of the synchronous time r. Although the trans- 
formation will change the explicit form of the coordi- 
nation bivector W" , this change involves only the time 
coordinate, and in the spatial averaging limit there is 
no difference between averages computed in the VPCs 
and those computed after the time redefinition. This 
again emphasizes the importance of the spatial averag- 
ing limit of spacetime averaging, if we are to succeed 
operationally in explicitly displaying the correlations as 
corrections to the standard cosmological equations. The 
correlation terms in Eqns. (|70p are therefore still inter- 
preted with respect to the volume preserving gauge, but 
are treated as functions of r. For the scale factor on the 
other hand, we have 



n 3 H- 1 da - H 
an — — — = _ni 

a dr 



FLRW 



-+3H 2 

a 



1 d 2 a 
a dr 2 



Further writing 

p(r) = p(i(r)) 
equations (p9"| become 
8nG N 



P (t) =p(t(r)), 



H 



FLRW 

lfa 
a dr 2 



-P- 7 



AttG 



N 



(p + 3p) + -S 



(i) 



(79) 



(80) 



(81a) 



(81b) 



We emphasize that the quantities Q^ 1 ** and defined 
in Eqns. (|70ap and (|70b[) as correlations in the volume 
preserving gauge, are to be thought of as functions of 
the synchronous time r, where the coordinate r itself 
was defined after the spatial averaging. Such an identi- 
fication is justified since we are dealing with scalar com- 
binations of these quantities. Note that and 
can be treated independently, apart from the constraints 
imposed by Eqn. (|41[) , which we turn to next. These 
conservation conditions can be decomposed into a scalar 
part and a 3- vector part, given respectively by 



v b Cl 



b;a 



o 



0. 



(82) 



In the synchronous gauge ((78)) for the FLRW metric, the 
scalar equation reads 



6H, 



FLRW 



G< x >) + (. 



d T SW +2H FLRW S^j =0. 

(83) 

We recall that this equation is a consequence of setting 
the correlation 3-form and the correlation 4-form to zero, 
and it relates the evolution of Q^ 1 ) and S^ 1 ' . The 3- vector 
equation (on imposing the first set of conditions in Eqn. 
([71])) simply gives d T C T A — 0, so that C T A = = constant, 
which also implies that C A = = constant and hence this 
equation gives nothing new. (We have used the relations 
C° = C T T , C° A = d 3 C T A and C A = {l/d 3 )C A where 
denotes the nonsynchronous time coordinate t.) 

The cosmological equations (f5Tj) . along with the con- 
straint equations (|77|) and (|83[) are the key results of this 
section. Subject to the acceptance of the volume pre- 
serving gauge on the underlying manifold M they can 
in principle be used to study the role of the correction 
terms resulting from spatial averaging. 



B. Results for an arbitrary gauge choice 

In this subsection, we will display the results obtained on 
assuming that the metric 



^ds 2 



-N 2 (t 7 x)dt 2 + h A B(t,x.)dx A dx 1 



(84) 



averages out to the FLRW metric in standard form with 
a nonsynchronous time coordinate t in general, to give 



^ds 2 



-f 2 (t)dt 2 + d 2 (t)S AB dx A dx I 



(85) 



In other words, we are assuming that the relations in 
Eqn. (|48| hold. Note that the averaging operation is no 
longer trivial, although we are still assuming an averag- 
ing on domains corresponding to "thin" time slices. We 
again split the averaged Einstein equations into scalar 
equations, and 3- vector and traceless 3-tensor equations. 
After transforming to the synchronous time coordinate 
r, now defined by 



f(t')dt' . 



(86) 



and again defining H = (l/d)(da/dt) and -Sflrw = 
(l/a)(da/dr) with a(r) = a(i(r)), the modified Fried- 
mann and Raychaudhuri equations read 



fi 2 

-"FLRW 



1 a a 
a dr 2 



8ttGn 



(87a) 



4ttG 



N 



(P- 



3p) + - 



75(i) +7 3(2) + £(2) 



(87b) 
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where the correlation terms are now defined using the 
relations, 



1 



r A pB 
1 OA 1 OB 



(88a) 



where the lower case index m in the last equation runs 
over all spacetime indices 0, 1,2,3, with the index re- 
ferring to the nonsynchronous time t. It is easy to check 
that Eqns. ([55]) reduce to Eqns. (fTTj) with the choice 
N = h~ x l 2 . The condition C b . a = has the scalar part, 



5(D 



>(2) 



(9 JK ) 



1 JB L KA I 



1 

T 2 

3H 
1 



r A F 



0A L 00 / 



2 (In/) 



V A V B \ 
1 JA L KB I 



H) 



JA l OK 



1 



5(2) _ _L (r A r°< 

° — p \ 00 1 AO 



(88d) 



We emphasize that averaging here refers to spatial av- 
eraging. Also (g JK ) = G JK = (l/a 2 )S JK , and the in- 
dex refers to the nonsynchronous time t. It is easy 
to check that and + V^ 1 correspond to cor- 
relations of (the bilocal extensions of) the extrinsic cur- 
vature with itself and with the time derivative of the 
lapse function. corresponds to correlations between 
the bilocal extensions of the Christoffcl symbols of the 
3-geometry, and S^ 2%> to correlations of the extension of 
the spatial derivative of the lapse function with itself and 
with the Christoffel symbols of the 3-geometry. Due to 
the way we have defined these correlations, one can also 
check that when the lapse function satisfies NVh = 1 (so 
that the averaging becomes trivial), we have "PC 1 ) = QW, 
5« = 5 (1) , and f> {2) = = S {2 \ where Q« and 5« 
were defined in Eqns. ([70)1 . The 3- vector and traceless 
3-tensor equations become 



1 

T- 



p0 pB 
1 OA 1 BO 



nU 1 AO ) 



r° r 

1 n r1 



+ ( 



9 JK ) 



p0 ?B 
1 ,/B 1 AK 



1 JA L BK I 



(89a) 



1 

T 2 



pA pB 
1 00 1 BO 



fBfA \ 
1 00 1 BO / 



+ <9 JK ) 



pi PB 
1 ,/B 1 OA' 



V A v 



b V 



./O- 1 BK J 

(89b) 







i r 

T 2 



v A r 

1 BO 1 0m / 



-pA pm \ 
1 mO 1 OB / 



7?(2) 



3 ds 



5(2)_5«-^(// 



3 ft (In/) 



(89c) 



(a T P (1) + 6i7 FLRW ^ (1) ) + (ftS (1) + 2ff FLR wS (1) 

+ 4ff FLRW (^ (2) +5 (2) ) =0, 

(90) 

while the 3- vector part, as before, gives nothing new and 
simply states d T C T A — 0. 

We can now state the main result of our paper in a clear 
and unambiguous manner, as follows : Having assumed 
that the FLRW spatial sections arise as the average of 
some gauge choice with lapse function N(t, x), spatial 3- 
metric fiAB{t, x) and shift vector N A set to zero for con- 
venience, we can construct the scalar quantities C^v v a 
and 7r^Cj + C b v b v a which, in coordinates natural to the 
FLRW metric take the form, 



L r p(D + 5(i) 



(91) 



C a b v b v a = - 



-K h a C a b + C%v b v a = 2 \pw + V (2) + S (2) l , 



with the various quantities being defined in Eqns. 
These scalars modify the usual cosmological equations as 
shown in Eqns. (J87J), and are themselves subject to the 
differential conditions (J5UJ- In addition, for consistency 
of our assumptions with the formalism, the underlying 
inhomogeneous metric is also subject to the conditions 



The combinations on the right hand sides of the re- 
lations (|9"Tj) can clearly be treated independently, apart 
from the conditions (|90|) . Further, since the correlation 



2-form has 40 independent components Z a bm l - n after im- 
posing all algebraic constraints, and since none of the 
four quantities V^- l \ V^ 2 \ and are trivially re- 
lated by these constraints, one can always treat these four 
functions independently of each other, subject only to 
the constraint in Eqn. (190|) . Before proceeding, we wish 
to make two remarks concerning the possible behaviour 
of the correction terms. It was mentioned in Ref. [l2j 
that assuming only spatial correlations, i.e. assuming 
that all components of the correlation 2-form Z a A with 
at least one index vanish, the corrections must be of 
the form of a spatial curvature term in the FLRW equa- 
tions. We can confirm this statement, since the above 
assumption amounts to setting and to 

zero, leaving only which must then evolve as ~ a~ 2 
because of Eqn. (|9H)) . Further, the main result of the 
Ref. [13] was that assuming the averaged metric to be 
FLRW in the conformal gauge (which corresponds to as- 
suming f(t) = a(t) in our case) and further assuming all 
components of the correlation 2-form to be constant, the 
corrections must again be of the form of a spatial curva- 
ture term. Eqns. (|88P show that this is indeed the case, 
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and Eqn. (|9Q[) shows that in this case one must have 
p(X) + + 5(2) _ 0, with no condition on the con- 
stant a 2 ^ 1 -* , which is consistent with a curvature term 
in the modified FLRW equations ([57)1 . Clearly though, 
the allowed behaviour of the correction terms which is 
consistent with Eqn. ([90]). is more general than that of a 
spatial curvature term, and it is not yet clear how these 
corrections behave in the real Universe. 

It is only fair to say that much of what we have de- 
scribed is already implicit in the work of Zalaletdinov 
and collaborators. What we have done here is to spell it 
out explicitly, emphasizing the relevance of spatial av- 
eraging for Cosmology. Further, the correction terms 
resulting from averaging have been displayed explicitly 
as scalars, which could be of help in applications and 
comparison with observations. A hitherto unappreciated 
fact which emerges is that even if some inhomogeneous 
geometry yields an FLRW Universe upon averaging over 
sufficiently large scales, this is not sufficient to guarantee 
consistency with the averaged Einstein equations. The 
consistency conditions (|59"|) must be satisfied - this fact 
could have potential significance in restricting the class 
of allowed initial perturbations in the early Universe, and 
should be investigated further. 

Another important goal of our paper is to attempt to 
compare the approaches of Zalaletdinov and Buchert, a 
topic to which we now turn. 



A comparison with the averaging formalism of 
Buchert 



The averaging formalism developed by Buchert is based 
exclusively on the manifold M. , and there is no analog of 
the averaged manifold M. in this scheme. Given an inho- 
mogeneous metric on M. one takes the trace of the Ein- 
stein equations in the inhomogeneous geometry, and car- 
ries out a spatial averaging of the inhomogeneous scalar 
equations. 

We recall in brief Buchert 's construction 0, by first 
writing down the averaged equations for the simplest 
case of pressureless and irrotational inhomogeneous dust. 
The metric can be written in synchronous and comoving 
gauge as 



ds 2 = -dt 2 + b AB (x,t)dx A dx 1 



(92) 



The Einstein equations can be split 3 into a set of scalar 
equations and a set of vector and traceless tensor equa- 
tions. The scalar equations are the Hamiltonian con- 
straint (|93aP and the evolution equation for 9 (|93b[) . 



-6 2 - 2a 2 = 16ttGp, 



K + d t O + e 2 = 12ttGp, 



(93a) 



(93b) 



where 1Z is the Ricci scalar of the 3-dimensional hyper- 
surface of constant t, O and a 2 are the expansion scalar 



and the shear scalar defined earlier and p is the inhomoge- 
neous matter density of the dust. Note that all quantities 
in Eqns. (f93|) generically depend on both position x and 
time t. Eqns. (|93a[) and (|93b[) can be combined to give 
Raychaudhuri's equation 



<9*e + \<d 2 + 2o- 2 + AirGp = . 



(94) 



The continuity equation dtp = —Op which gives the evo- 
lution of p, is consistent with Eqns. (I93aj) . (I93bj) . Only 
scalar Einstein equations are considered, since the spa- 
tial average of a scalar quantity can be defined in a gauge 
covariant manner, within a given foliation of space-time. 
We return to this point below. For the space-time de- 
scribed by (f9"2")l . the spatial average of a scalar ^(x, i) 
over a comoving domain T> at time t is defined by 



v 



1 

Vv 



d 3 xVb 1 $> . 



(95) 



where b is the determinant of the 3-metric dab and 
Vv is the volume of the comoving domain given by 
Vd = J v d 3 xVb. Spatial averaging is, by definition, not 
generally covariant. Thus the choice of foliation is rele- 
vant, and should be motivated on physical grounds. In 
the context of cosmology, averaging over freely- falling ob- 
servers is a natural choice, especially when one intends 
to compare the results with standard FLRW cosmology. 
Following the definition (|95[) the following commutation 
relation then holds [3j 

dt - (d t *) v = (*e) B - (e) v , (96) 

which yields for the expansion scalar 

d t (G) v - (d t Q) v = (Q 2 ) v - (0) 2 V . (97) 
Introducing the dimcnsionless scale factor a-p = 

1 /3 

(Vu/Vui) normalized by the volume of the domain 
T> at some initial time t,, we can average the scalar Ein- 
stein equations (|93a[) . (|93b[) and the continuity equation 
to obtain [3[ 



d t ( P ) v = - (e) v ( P ) 



V > 



(G) v = 3 



d t av 



(98a) 



d t a v 



up 



8ttG 



G 



(Qv + (K) v ) , (98b) 



dfa-v 



^(p)v + \Qv- (98c) 



Here, the 'kinematical backreaction' Qx> is given by 

Qv^l((0 2 ) v -{Q) 2 V ) -2 {a 2 ) v (99) 
and is a spatial constant over the domain T>. 
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A necessary condition for (|98c|) to integrate to (|98b|) 
takes the form of the following differential equation in- 
volving Qx> and (TZ) V , 

d t Qv + 6—Qv + d t (Jl) v + 2^(K) v = 0. (100) 
a-p av 

The equations above describe the essence of Buchert's 
averaging formalism, for the dust case. We note that 
the remaining eight Einstein equations for the inhomo- 
geneous geometry, which are not scalar equations, are not 
averaged. These are the five evolution equations for the 
trace- free part of the shear, 

d t (a£) = -eai -Ki + -^i (a 2 - ^6 2 + SnGp^j . 

(101) 

and the three equations relating the spatial variation of 
the shear and the expansion, 

"B||A = f e ll*- ( 102 ) 

Here, lZg is the spatial Ricci tensor and, in Buchert's 
notation, a 1 1 denotes covariant derivative with respect 
to the 3-metric. 

In analogy with the dust case, Buchert's averaging for- 
malism can be applied to the case of a perfect fluid 
by starting from the metric 

ds 2 = -N 2 dt 2 + b AB dx A dx B . (103) 

The averaged scalar Einstein equations for the scale fac- 
tor ax> are 

3 ^^ + 4 7r G(jV 2 0) + 3 P )> 1J = S 7 , + P 7 ,, (104) 



6Hl-16wG(N 2 p) v = -Q V -(N 2 K) V ; H v = , 

(105) 

where the kinematical backreaction Qt> is given by 

Qv = \ (((m 2 ) v - 2 {N 2 a 2 ) v , (106) 

and the dynamical backreaction Vx> is given by 

Vv = (N 2 A) V + (@d t N) v , (107) 

where A = V j(u*V 'jU J ) is the 4-divergence of the 4- 
acceleration of the fluid. Eqn. (|105[) follows as an integral 
from Eqn. (|104ll if and only if the relation 

d t Q v + MvQv + d t {N 2 K) V + 2H V (N 2 K) V + AH V V V 

- 16^G [d t (N 2 p) v + 3H V (N 2 (p + p)) v ] = , 

(108) 

is satisfied. There are also the unaveraged equations 
(which we do not display here) for the shear, analogous 
to the shear equations (|101|) and (|102p for dust. 



Buchert's approach is the only other approach, apart 
from Zalaletdinov's MG, which is capable of treating in- 
homogeneities in a nonperturbative manner, although it 
is limited to using only scalar quantities within a chosen 
3+1 splitting of spacetime. Buchert takes the trace of the 
Einstein equations in the inhomogeneous geometry, and 
averages these inhomogeneous scalar equations. In the 
context of Zalaletdinov's MG however, we have used the 
existence of the vector field v a in the FLRW spacetime 
to construct scalar equations after averaging the full Ein- 
stein equations. As far as observations are concerned, it 
has been noted by Buchert and Carfora [13] , that the spa- 
tially averaged matter density {p) v defined by Buchert 
is not the appropriate observationally relevant quantity 
- the "observed" matter density (and pressure) is actu- 
ally defined in a homogeneous space. Since we have done 
precisely this in Eqn. (|45p . we are directly dealing with 
the appropriate observationally relevant quantity in the 
MG framework. 

Another important difference between the two ap- 
proaches is the averaging operation itself. Buchert's spa- 
tial average, defined for scalar quantities, is given (for 
some scalar ^(t, x A )) by ([55)) above. On the other hand 
the averaging operation we have been using (given by 
Eqn. (|53[) using the volume preserving gauge) is a limit 
of a spacetime averaging defined using the coordination 
bivector W° , and is different from the one in Eqn. (|95[) . 

Most importantly though, Buchert's averaging scheme 
by itself does not incorporate the concept of an averaged 
manifold M. (although the work of Buchert and Carfora 
[I?! does deal with 3-spaces of constant curvature). In 
a recent paper j2Q| we had argued that Buchert's "effec- 
tive scale factor" ax>{t) = (Vv(t) /Voi^m)) 1 ^ must be the 
scale factor for the metric of the averaged manifold, upto 
some corrections arising due to such effects as calculated 
by Buchert and Carfora. In the present work however, 
it is clear that such a suggestion is necessarily incom- 
plete due to the presence of Eqns. (|77|) constraining the 
underlying geometry. These constraints are in general 
nontrivial and hence indicate that it is not sufficient to 
assume that the metric of the inhomogeneous manifold 
averages out to the FLRW form - there are additional 
conditions which the correlations must satisfy. 

To our understanding, Buchert's averaging formalism 
is a valid aproach, even though it is based on a spatial 
averaging. A central difference from the MG approach 
is the issue of closure : not all the Einstein equations 
have been averaged in Buchert's approach, but only the 
scalar ones. This puts a constraint on the allowed solu- 
tions considered for the averaged equations: (|98p for the 
dust case, and (|104|) and (|105|) for the fluid case. So- 
lutions to these equations must necessarily be checked 
for consistency with the unaveraged equations for the 
shear. Further, averaging over successively larger scales 
can bring in additional corrections to the averaged equa- 
tions, as discussed by Buchert and Carfora. Also, if one 
does not wish to identify Buchert's a-p with the scale fac- 
tor in FLRW cosmology, one is compelled to develop a 
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whole new set of ideas in order to try and compare theory 
with observation. On the other hand, if one does iden- 
tify a-p with the scale factor, comparison with standard 
cosmology becomes more convenient, but this brings in 
additional constraints on the underlying inhomogeneous 
geometry. Thus our conclusion is that the Buchert for- 
malism is a correct and tractable averaging scheme, pro- 
vided all the caveats pointed out in this paragraph are 
taken care of. Also, when these caveats have been taken 
care of correctly, the Buchert formalism is expected to 
give the same physical results as the MG approach. We 
recall that in the covariant MG approach also, once a 
spacetime geometry has been identified for the averaged 
manifold M, a gauge must be selected for the geometry 
on the underlying manifold, in order to explicitly com- 
pute the correction scalars for comparison with observa- 
tion. 

The advantage of the MG approach is that it accom- 
plishes in a neat package what the Buchert approach, 
with its attendant caveats, sets out to do. In the MG 
approach, there are no unaveraged shear equations, be- 
cause the trace of the Einstein equations has been taken 
after performing the averaging on the underlying geom- 
etry. Since the averaged geometry is FLRW, the shear 
is zero by definition. There is a natural metric on the 
averaged manifold by construction, the FLRW metric. 
The correlations satisfy additional constraints, given by 
Eqns. (|77| . Thus, once a gauge has been chosen and if 
one can overcome the computational complexity of the 
averaging operation, the cosmological equations derived 
by us in the MG approach are complete and ready for 
application, without any further caveats. 

In spite of these differences, our equations (f8Tj) and 
([83]) for the volume preserving gauge are strikingly sim- 
ilar to Buchert 's effective FLRW equations and their in- 
tegrability condition in the dust case; and in the case 
of general N, the role of Buchert's dynamical backreac- 
tion Vv in Eqns. p04|) and (|108j) is identical to that 
of our combination of ("P (2) + 5 (2) ) in Eqns. ((87b)) and 
(|90f . Concentrating on the volume preserving case, the 
structure of the correlation is identical to Buchert's 
kinematical backreaction Qt> (or Qx> in the general case). 
The correlation <SW appears in place of the averaged 3- 
Ricci scalar (1Z) V in Buchert's dust equations. This is 
not unreasonable since Buchert's (72) -p can be thought 
of as (72.) p = &kx>/a\ ) -\- corrections, where Qk-o/a^, rep- 
resents the 3-Ricci scalar on the averaged manifold which 
in our case is zero, and hence S 1 ^ represents the correc- 
tions due to averaging. Further, these similarities are in 
spite of the fact that our correlations were defined as- 
suming that a volume preserving gauge averages out to 
the FLRW 3-metric in standard form, whereas Buchert's 
averaging is most naturally adapted to beginning with 
a synchronous gauge. This remarkable feature, at least 
to our understanding, does not seem to have any deeper 
meaning - it simply seems to arise from the structure 
of the Einstein equations themselves, together with our 
assumption D^Z"^- = 0. In the absence of this lat- 



ter condition, one would have to consider the correlation 
3- and 4-forms mentioned earlier, and the structure of 
the correlation terms and their "conservation" equations 
would be far more complicated. 

An entirely different outlook towards his approach has 
been emphasized to us by Buchert [2^ |. According to 
Buchert, the absence of an averaged manifold M. is not 
to be thought of as a 'caveat', but as a feature deliber- 
ately retained 'on purpose'. The actual inhomogeneous 
Universe is regarded by Buchert as the only fundamental 
entity, and the introduction of an averaged Universe is in 
fact regarded as an unphysical and unnecessary approx- 
imation. As we mentioned earlier, this is probably the 
most important difference between MG and Buchert's ap- 
proach. In the latter, contact with observations is to be 
made by constructing averaged quantities, such as the 
scalars defined earlier in this section, and by introduc- 
ing the expansion factor ap. The assertion here is that 
the averaging of geometry, as discussed in MG or in the 
Renormalization Group approach of Buchert and Car- 
fora [17] ] is not an indispensable step in comparing the 
inhomogeneous Universe with actual observations. The 
need for averaging of geometry is to be physically sep- 
arated from simply looking at effective properties (such 
as the constructed scalars) which can be defined for any 
inhomogeneous metric. Averaging of geometry becomes 
relevant if (i) an observer insists on interpreting the data 
in a FLRW template model, so that (s)he needs a map- 
ping from the actual inhomogeneous slice and its average 
properties to the corresponding properties in this tem- 
plate, or (ii) one desires a mock metric, to sort of have a 
thermodynamic effective metric to approximate the real 
one. (In this context it should perhaps also be mentioned 
that the importance of a thin time-slice approximation of 
spacetime averaging (as opposed to a strict spatial aver- 
aging) has been stressed also by Buchert [|[.) 



V. DISCUSSION 

In this paper we have addressed the issue of modifica- 
tions to the standard cosmological equations arising out 
of explicitly accounting for the averaging procedure that 
must necessarily be performed when studying Cosmology. 
While this issue has been dealt with in the literature in 
various forms, the effect of such modifications is still a 
subject of debate, mainly due to ambiguities in the aver- 
aging schemes available. We have applied the formalism 
of Zalaletdinov's Macroscopic Gravity (MG) which is a 
fully covariant and nonperturbative averaging scheme, in 
an attempt to construct gauge independent corrections 
to the standard FLRW equations. 

We find that one cannot escape the problem of gauge 
dependence entirely, which is mainly due to the fact that 
Weyl's postulate which is commonly used in the standard 
approach to Cosmology, must be refined in the context 
of averaging in General Relativity, and one is forced to 
assume that the FLRW spatial sections in their natu- 
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ral coordinates arise from averaging the inhomogeneous 
metric in a particular gauge. However, the choice of this 
gauge can at least formally be left unspecified, and space- 
time scalar corrections to the FLRW equations can be 
constructed. This partially removes the criticism usually 
faced by adherents of the averaging approach, which es- 
sentially states that effects of averaging are likely to be 
gauge artifacts. We have shown that these effects are 
independent of the choice of coordinates in the averaged 
manifold, and depend on one gauge choice - a choice 
which itself is fundamental to the assumption of large 
scale homogeneity and isotropy. 

One issue which we did not address was the fact that 
the scale of averaging is likely to be different at differ- 
ent epochs. The scale Lflrw ~ 100/i _1 Mpc which we 
mentioned in Section IIII1 would be appropriate for the 
present day Universe. For the early Universe however, 
this scale will probably be very different. Nevertheless, 
since the formalism simply assumes the existence of a 
scale at which homogeneity sets in, without being af- 
fected by its actual value, the equations we have derived 
within the MG framework will apply to both the present 
day as well as to the early Universe, albeit with different 
averaging scales. If one assumes that the length scale of 
averaging varies slowly compared to the Hubble expan- 
sion scale, then it would be appropriate to apply the for- 
malism independently to say, Supernovae data and the 
last scattering epoch. However, this issue needs to be 
dealt with more precisely than discussed here. 

This brings us to the most important question to be 
addressed : How does one apply this formalism without 
knowing the "true" inhomogeneous metric of the Uni- 
verse? One possible line of attack may be as follows : 
In the early Universe at least, it might be reasonable to 
assume that the "true" metric is a perturbation around 
the FLRW metric. Using some standard gauge, making 
the necessary assumptions and ensuring that all required 
conditions are satisfied, one should then be able to ex- 
plicitly compute the scalar corrections given the initial 
power spectrum of the fluctuations. The behaviour of 
these corrections (in the perturbative regime) may yield 
some insight as to how the corrections might evolve in 
the present epoch (nonlinear regime). One could then 
construct reasonable models for the functions S^ 1 ', 
V {2) and 5 (2) defined in Section HVBl Note that in this 
regime, one needn't worry about the explicit construc- 
tion of these quantities - what one needs are physically 
reasonable models of the time evolution of these quanti- 
ties (see, e.g., Ref. [Hj]). We hope to commence such a 
program in the near future. 
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APPENDIX A 

In this appendix we give proofs of several results that 
were used in the text. 



1. Analysis of T>n9 = 

We start with the metric 

( M Us 2 = g o(t,x)dt 2 + g AB (t,x)dx A dx B , (Al) 

on A4 and assume that it averages out to the FLRW form 
(Eqn. O): 

Goo = (3oo) = -f(t) ■ G OA = (g OA )=0; 
Gab = (sab) = a 2 (t)S AB . (A2) 



We will analyze the second relation of Eqn. (|31j) and 
show that it leads to the result £/ y = g 13 — G % i — 0, where 
tt a b refers to the connection 1-forms associated with Gy . 
For this section we use the notation H = (l/a)(da/dt). 
We have 



dg ab + n a j gi b + n b ,g aj = . 



(A3) 



Consider the three cases (a = b = 0), (a = 0, b = B) and 
(a = A, b — B) in turn. From the first case (a = b = 0) 
we can conclude that 



3 UU (*,x) 



fc(x) 



d A k(x) = 2a 2 H6 AB g 



OB 



(A4a) 



(A4b) 



where fc(x) is a positive definite function (so that the 
metric signature is preserved) which arises as an integra- 
tion constant and is constrained by Eqn. (|A4bp . The 
second case (a = 0, b = B) leads to 



-OB 



m B (x) 



i-3W(x) + ^HS AJ g AB - — 



(A5a) 



HSj = . (A5b) 



where m s (x) is a 3- vector that arises as a constant of 
integration like fc(x), and is constrained by Eqn. (|A5b[) . 
Finally, the last case (a = A, b = B) leads to 



-AB 



a 2 (t) 



S ^(x). 



(A6a) 
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^dj S AB (x) + jjH (<5> s (x) + 5j m A (x)) 



0. 



(A6b) 

Here s AB (x) is another constant of integration, a sym- 
metric 3-tensor. It can be easily argued that since / and 
a are not a priori related, both sides of Eqn. (IA4b[) must 
vanish, which immediately tells us that the vector m B (x) 
must vanish, and the function fc(x) must be a constant, 

fc(x) = k= constant ; m s (x) = . (A7) 

Equations (|A5b|) and (|A6b[) then give us 

s AB (x) = kS AB , (A8) 

with the same constant k as in Eqn. (|A7|) . Finally, 
putting everything together we find 



=00 



k 



9° A = 



-AB 



& X AB 



kG K 



(A9) 



The constant k is not constrained by any of the equa- 
tions and appears to be a free parameter in the theory. 
The modified Einstein equations (|35[) show that k can 
be absorbed into the averaged energy momentum ten- 
sor. We will for simplicity assume k to be unity thereby 
obtaining, as required 

U lj = g lj - G lj = . (A10) 

2. Analysis of the condition (r£ c ) = ( FLHW >r£ c 

Here we will assume that the line element on M. is in the 
volume preserving gauge 

dP 



A,R 



+ h AB (t,x)dx dx 



(All) 



so that the averaging is trivial, and the metric and aver- 
ages out to the FLRW line element on M. given by 



(^W = - 



dp 



\t)5 AB dx A dx B 



(h)(t) 



(A12) 



where we used the condition (l/h) = 1/ (h) that follows 
from g 00 = G 00 . The conditions (rg c ) = ( FL R w )r° c then 
result in the following : 



r° 

1 oo ■ 




(A13a) 


1 OA ■ 


(d A ()jxVh))=0, 


(A13b) 


r A . 

1 oo ■ 


1 h AB \ 
/—d B (lnVh)j=0, 


(A13c) 


r A . 
1 0B • 




(A13d) 


r o 

1 AB ■ 


(yhQA B ) = (h) a 2 H6 AB , 


(A13e) 


r A 

1 BC ■ 




(A13f) 



Eqns. (IAl3bj) and jAilf]) are consistent with each other 
since {3) T^ A = d B (lnVK), and Eqn. (fAT3c|) is consis- 
tent with the assumption Eqn. (|30a|l . The trace of Eqn. 

(|A13d|) gives ((l/Vh)Q) = 3if. However, using Eqn. 

([rJTjl we have (l/ Vh)Q = NO = dt{lnVh), and com- 
bined with Eqn. (|A13a,p this gives 



-a f (ln(/ l ))=39 f (lna) 



(h) = a 6 



(A14) 



where we have set an arbitrary proportionality constant 
(representing rescaling of the time coordinate by a con- 
stant) to unity. This establishes the last equality in Eqn. 



Finally, consider the trace ((h AB ) / (h)) (^fhQ AB 

using the condition g AB = G AB , Eqn. (|A13e|) and the 
trace of Eqn. (|Al3dj) , this gives us 



(h AB ) I r 



AB 



SjAB , 

^(VhO AB 



(h) 

=3H=( 7n e) = 



-.AB 



-{VhQ AB ) 
(A15) 



On using the condition (|30a|) this leads to 



' (h AB ) /h 
(h) 



,AB ' 



(r%) = o, 



(A16) 



which is consistent with the assumption 



(h AB ) _/h 
(h) 



,AB 



(A17) 
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